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摘 要 : 在 实际 通信 环境 中 ， 由 于 传播 环境 的 复杂 性 使 空间 中 存在 大 量 的 相干 信号 ， 从 而 导致 信 源 协 方差 矩阵 的 秩 
亏 缺 。 为 使 得 矩阵 的 秩 恢复 到 等 于 信号 源 数 并 解决 相干 信 源 波 达 方向 (direction of arrival ，DOA) 估计 问题 ， 提 出 
了 一 种 混合 型 MUSIC 算法 。 po 处 理 ， 并 将 得 到 的 新 协 方差 志 本 起 
阵 应 用 于 改进 的 IMUSIC 算法 进行 信号 数据 处 理 分 析 ， 得 到 相干 信号 的 DOA 角度 估计 。 仿 真 结果 表明 ， 在 信 嗓 比 
低 的 情况 下 ， 信 和 号 We i i 能 准确 地 估计 出 信 源 的 DOA， 验 证 了 该 算法 的 
高 分 辨认 和 高 性 能 。 

关键 词 : 相干 信 源 ; 多 重信 号 分 类 ; 空间 平滑 ; 波 达 方向 估计 
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Doa estimation of coherent sources based on hybrid music algorithm 


= Zhang Quanl, Li Simin?, Tang Zhiling'!i 

by (1. Institute of Electrical Engineering & Automation, Guilin University of Electronic Technology, Guilin Guangxi 541004, 
China; 2. Institute of Electrical & Information Engineering, Guangxi University of Science & Technology, Liuzhou Guangxi 
545006, China) 


Abstract: In the practical communication environment, because of the complexity of the propagation environment, there are 
a lot of coherent signals in the space, which leads to the rank defect of the source covariance matrix. In order to restore the 
rank of the matrix to be equal to the number of signal sources and solve the problem of DOA estimation, this paper proposed 
a hybrid MUSIC algorithm. It pre-processed the antenna array by the forward and backward spatial smoothing technique, 
and applied the new covariance vector matrix to the improved IMUSIC algorithm for signal processing and analysis, and 


obtained the DOA angle estimation of coherent signals. The simulation results show that the hybrid MUSIC algorithm can 
accurately estimate the DOA of the source when the signal-to-noise ratio is low, the signal interval is small and the 
correlation signal exists, which verifies the high resolution and high performance of the algorithm. 


Key words: coherent signal; multiple signal classification; spatial smoothing; direction of arrival estimation 


GG 0 引言 DOA 估计 Root-MUSIC 算法 , 它 将 空间 谱 峰 值 搜索 过 程 转换 
四 全 为 多 项 式 求解 过 程 ， 该 算法 测 向 精度 高 ， 处 理 小 样本 性 能 优 
CS 空间 信和 号 波 达 方 向 估计 是 阵列 信号 处 理 中 的 关键 问题 ， 越 。 文 献 [6] 提 出 改进 MUSIC 算法 将 信号 协 方差 矩阵 进行 预 
并 用 于 雷达 、 声 呐 以 及 无 线 通信 等 许多 场合 ， 可 以 实现 对 目 处 理 ， 使 信号 协 方差 矩阵 分 解 得 到 的 噪声 子 空间 与 信号 子 空 
标的 定位 与 跟踪 和 信号 的 空间 滤波 。 闻 正 交 ， 降 低 了 噪声 干扰 。 文 献 [7] 提 出 使 用 最 大 似 然 方法 来 
在 众多 国内 外 学 者 的 研究 探索 下 ， 不 同 的 DOA 估计 算 推导 阵列 的 加 权 MUSIC (w-MUSIC) 算法 ， 其 通过 组 合子 
法 相继 被 提出 。 文 献 [1] 首 次 提出 多 信号 分 类 (multiple signal 。” 阵列 的 加 权 MUSIC 谱 来 获得 整体 空间 谱 。 文 献 [8] 提 出 基于 
classification，MUSIC ) 算法 ， wo 以 信号 方向 向 量 与 噪声 特征 空间 的 DOA 估计 算法 ， 该 算法 结合 空间 平滑 技术 并 利 
子 空间 的 正 交 特性 为 基础 ， 构 造 空间 谱 函 数 通 过 谱 峰 搜索 进 用 阵列 协 方差 矩阵 的 一 个 子 矩 阵 快速 得 到 信号 子 空间 。 无 须 
行 DOA 估计 。MUSIC 算法 在 4 下 有 具有 很 高 的 分 辨 ”特征 值 分 解 或 奇异 值 分 解 , 降低 算法 的 复杂 度 并 提高 了 DOA 


ml 


力 、 估 计 精 度 及 稳定 性 ， 并 适用 于 各 种 形式 的 阵列 天 线 。 之 十 计 效 率 。 

后 在 MUSIC 算法 的 基础 上 , 演变 出 基于 特征 空间 的 DOA 估 针对 相干 信号 会 导致 信 源 协 方差 矩阵 的 秩 亏 缺 ， 使 得 常 
计算 法 、 修 正 MUSIC 算法 、 改 进 MUSIC 算法 、 求 根 MUSIC 用 的 信号 子 空间 类 算法 在 DOA 估计 过 程 中 失效 。 目 前 解 光 
算法 (Root-MUSIC)〉 和 加 权 MUSIC 算法 等 。 秩 亏 缺 的 问题 ， 多 数 方法 是 使 信号 协 方差 矩阵 的 秩 得 到 有 效 

文献 [2] 提 出 修正 MUSIC 算法 ， 其 实质 是 取 前 后 向 空间 恢复， 或 重新 构造 一 个 秩 等 于 信号 源 数 的 矩阵 。 

平滑 算法 中 子 阵 个 数 为 1 的 情况 ， 使 得 子 阵 长 度 等 于 阵 元 个 空间 平滑 算法 由 Evans 等 人 [2 首先 提出 ， 其 核心 原理 
数 。 该 算法 不 会 减 小 阵列 的 孔径 ， 在 信 品 比较 低 的 情况 下 依 、 ”是 将 等 距 线 阵 分 成 若干 个 相互 重 琶 的 子 阵列 ， 各 子 阵列 的 阵 
| 有 较 高 的 分 辩 力 。 文 献 [3~5] 提 出 一 种 基于 特征 分 解 的 ” 列 流 形 相同 ， 而 各 子 阵列 的 协 方差 矩阵 可 以 进行 平均 运算 ， 
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其 成 为 有 效 预 处 理 相 干 信 源 的 算法 。 文 献 [11] 使 用 多 个 循环 
相关 MCC) 和 空间 平滑 处 理 的 同 相 组 合 对 信号 进行 DOA 
估计 。 该 方法 在 使 用 两 个 循环 相关 和 矩阵 和 相应 的 循环 后 向 矩 
阵 (2MCC-CB) 情况 下 估计 精度 最 高 ， 增 强 了 DOA 估计 性 
能 。 文 献 [12] 提 出 前 后 向 空间 平滑 (FBSS) 算法 ， 首 先 将 前 
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实现 去 相干 。 文 献 [11~13] 通 过 对 空间 平滑 算法 进行 改进 ， 使 


为 0、 方 差 为 om 的 白 噪声 , 且 与 信号 源 不 相关 。 
4=[a(9),a(@)…,a(9,)] 是 方向 矩阵 。 其 中 a(4) 是 均匀 线 
阵 响应 矢量 可 表示 为 


了 
a(0)= [ exp[ -2rdsng - ep| -Pr -0D2sin 可 


并 


避 平 滑 的 子 阵列 共 思 倒置 后 形成 后 向 空间 平滑 ， 再 与 前 向 平 
滑 结 合 形成 前 后 向 空间 平滑 算法 。 空 间 平滑 技术 及 其 改进 算 
法 对 相干 信号 源 的 DOA 估计 具有 理想 的 性 能 ， 且 计算 量 很 
小 ， 便 于 实现 。 但 它 是 通过 牺牲 有 效 阵 元 数 来 换取 的 ， 对 阵 
列 孔 径 有 一 定 的 损失 ， 且 在 低 信 
[14~15] 提 出 了 一 种 改进 的 前 向 /后 向 空间 平滑 Root-MUSIC 
算法 ， 通 过 利用 了 实 值 协 方差 矩阵 的 特征 值 分 解 (EVD)， 
降低 了 Root-MUSIC 算法 的 运算 复杂 度 。 文 献 [16] 提 出 一 种 
改进 的 加 权 空 间 平 滑 (CIWSS ) 技 术 用 于 相干 信 源 DOA 估计 。 
该 方法 通过 应 用 舱 套 空间 平滑 算法 来 获得 权重 矩阵 , 与 WSS 
方法 相 比 在 没有 先 验 信息 和 初步 估计 角度 值 的 情况 下 也 能 获 
得 较 高 的 分 辨 率 和 估计 性 能 。 文 献 [17] 利 用 FB 技术 提出 了 一 
种 实 值 MUSIC 算法 ， 将 测 向 过 程 转换 为 真实 域 ， 从 而 显 着 
降低 计算 量 。 文 献 [18] 提 出 了 LMMUSIC 算法 ， 将 前 后 向 空 
间 平 滑 算法 与 MMUISC 算法 的 结合 并 应 用 于 车 载 防 撞 雷 达 。 
鉴于 以 上 分 析 ， 本 文 在 经 典 MUSIC 算法 和 当前 提出 的 
解 相 干 的 算法 基础 上 ， 提 出 混合 型 MUSIC 算法 。 该 算法 将 
前 后 向 空间 平滑 算法 和 改进 的 IMUSIC 算法 进行 融合 ， 并 通 
过 MATLAB 仿真 将 其 与 前 后 向 空间 平滑 和 改进 的 IMUSIC 
算法 以 及 传统 MUSIC 算法 进行 对 比 。 仿 真 结 果 表 明 ， 混 合 
型 MUSIC 算法 在 处 理 相干 信 源 DOA 估计 时 具有 较 高 的 测 向 
精度 和 分 辨 率 。 
1 ”问题 描述 
1.1 均匀 线 阵 信号 模型 

如 图 1 所 示 ， 由 M 根 天 线 阵 元 均匀 直线 排列 ， 单元 间距 
为 a ， 信 号 的 波长 为 4， 布置 成 一 个 阵列 天 线 。 假 设 有 
P(P<M) 个 窄带 信号 源 平 面 波 辐射 到 线 阵 上 ， 信 源 波 达 方 向 
8(i=42…,p) ， 则 第 :次 采样 点 数 的 信号 模型 可 表示 为 

X(t)=AS(1)+N (7) (1) 

其 中 : 于 (=[(0 (zw (1)] 表示 M 个 天 线 单元 数据 矩阵 ; 
7 表示 转 


SO)=[s(),s(0)…sw(O] 是 信 源 矩阵 ,其 中 8(D) 为 第 i 个 平面 
波 的 复 振幅 。 


下 


®@' © © © 

M 4 3 2 
所 

图 1 均匀 天 线 线 阵 示 意图 


Fig. 1 Schematic diagram of uniform antenna line array 


NOD)= 思 (0)( 人 mx( 是 噪声 矩阵 ,其 中 m(D) 是 均值 


其 中 4 是 信 源 i 的 波 达 方向 ,:=12…,P。 
1.2 MUSIC 算法 原理 

根据 上 述 均 匀 线 阵 信 号 模型 ， 由 式 (1) 可 求 得 阵列 信和 号 
的 协 方差 矩阵 为 
R=E[X(i)X"(1)]=AE[S()S" (A +o 
R=ARsA” + ol 2 
其 中 : Rs 是 信号 的 协 方差 矩阵 ;c27w 为 噪声 协 方差 矩阵 ;六 
为 W 阶级 单位 矩阵 。 对 RK 进行 特征 分 解 有 
R=UsEs UY +UVENUN (3) 
其 中 : Us 是 由 大 特征 值 对 应 的 特征 矢量 张 成 的 信号 子 空间 ; 
Us 是 由 小 特征 值 对 应 的 特征 矢量 张 成 的 噪声 子 空间 ; Z; 为 
?个 大 特征 值 构 成 的 对 角 阵 。 
根据 方向 矩阵 4 中 的 各 个 列 向 量 与 噪声 子 空间 正 交 ， 故 


有 
Uga(0)=0 (4) 
考虑 到 实际 接收 数据 矩阵 是 有 限 长 的 ， 即 数据 协 方差 矩 

阵 的 最 大 似 然 估计 为 


a。 1 之 本 
Rs (5) 
对 月 进 行 特征 分 解 可 得 到 噪声 子 空间 特征 矢量 矩阵 ， 
噪声 特征 向 量 和 信号 向 量 的 正 交 关系 , 得 到 阵列 空间 谱 函数 : 


1 
Puc (0) = (DORa 00) 


式 (6)， 当 a(9) 与 ,的 各 列 正 交 时 分 母 为 零 , 但 由 于 


噪声 的 存在 ， 使 得 Rusc 存在 尖峰 ， 通 过 寻找 峰值 来 估计 到 达 
角 。 
1.3 ”相干 信 源 对 DOA 估计 影响 

下 面 根据 相干 信号 源 数学 模型 ， 研 究 相 干 信 源 对 DOA 
估计 影响 。 对 于 两 个 平稳 信号 $:(?) 和 5 (0 ,定义 它们 的 相关 


OO) 

"alls oF efls, oF] 07) 
施 瓦 兹 不 等 式 可 知 ， 相 关系 数 满足 1Py 1 。 因 此 ， 对 
于 不 同 信号 之 间 的 相关 性 可 表示 为 当 Pp)=0 时 ， 称 5,() 和 
SO) 不 相关 ; 当 0d4mkl 时 ， 称 SC) 和 Si(D 部 分 相关 ， 当 
|py 上 = 时 ， 称 S;(1) 和 Sj (7) 相干 o 
信号 源 相干 时 其 数学 表示 为 :， 相 
常数 。 假 设 设 有 个 相干 信 源 ， 即 
Si(1)= S50(1), i=1,2,.…,n (8) 
式 中 ;Su(D 为 生成 信 源 ， 因 为 它 生成 了 入 射 到 阵列 上 的 # 个 
相干 信号 源 。 故 可 得 相干 信号 源 数学 模型 : 


Ql 


I 


t 


I 


信 源 间 只 差 一 个 复 


X(D=4| 一 |so(D)+N(GD)=4pSu(D+NO) (9) 
Cn 
其 中 : 2 为 复 常 数组 成 的 ax1 维 矢量 。 根 据 相干 信 号 源 模型 ， 
假设 有 个 入 射 信号 源 S17),5,(7)…,5, (7?) , 则 阵列 信号 可 表示 
为 
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SO 
KU)=[c(g) a(@) … $200 


a(0)] ”|+NO) 


(10) 
S, (1) 
信 源 自 相关 和 矩阵 为 R=ELXX”]=A4Rs;4”+o* ， 则 信号 源 
相关 矩阵 为 


SODSTCD) S(t)S3(1) ~ 
pS2(0S82() SS 


Si(1)S, (1) 
S, (1) 5, (7) 
S, (0)S:() S$, (ns: (1) 5 (Ds) 
of DO 7 Puoio, 


P2nO2On 


(11) 


paion 03 


# 2 
DCIO， DOazO， 人 On 


5S,(1) 的 功率 。 当 个 信 


其 中 :oo 为 信号 $1(?),S,(?),…， 


号 完全 相干 时 ， 其 中 i 和 j 都 小 于 等 于 n。 阵 列 接收 


的 数据 协 方差 R; 的 秩 降 为 1， 导致 信号 子 空间 的 维 数 小 于 信 
号 源 数 ， 从 而 无 法 正确 估计 信号 源 方向 。 因 此 针对 含有 相干 
言 号 源 的 DOA 估计 核心 问题 是 将 信号 协 方差 矩阵 Rs 的 秩 恢 
复 到 等 于 信号 源 的 个 数 P 。 
2 ”改进 算法 

混合 型 MUSIC 算法 实质 是 前 后 向 空间 平滑 算法 和 改进 
的 IMUISC 算法 的 融合 。 首 先 通过 对 天 线 中 的 阵 元 分 别 进行 
前 向 和 后 向 平滑 预 处 理 后 ， 使 得 阵列 协 方差 矩阵 的 秩 恢复 到 
言 号 源 个 数 ， 将 得 到 的 新 协 方差 矢量 矩阵 按照 改进 的 
IMUSIC 算法 进行 信号 数据 处 理 分 析 ， 得 到 相干 信号 的 DOA 
角度 估计 。 下 面 对 该 算法 的 阵列 预 处 理 以 及 具体 实现 步 又 进 
行 说 明 。 
前 后 向 空间 平滑 算法 把 前 向 空间 平滑 与 共 斩 后 向 空间 


lp;l=1 9 


以 
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滑 算 法 结合 ， er 的 。 阵 列 预 处 理 
部 分 是 通过 前 后 向 空间 平滑 算法 实现 的 。 两 个 相干 信号 同时 
进入 天 线 阵列 时 ， 会 使 得 阵列 协 方差 矩阵 的 信号 子 空间 的 矢 
量 减 少 : 但 若 两 个 相干 信号 同时 进入 不 同 的 天 线 阵 列 ， 那 么 
这 两 个 阵列 协 方差 矩阵 之 和 的 信号 子 空间 的 矢量 数 等 于 目标 
数 。 基 于 这 种 思想 ， 利 用 均匀 线 阵 的 平移 不 变形 ， 将 等 距 线 
阵 划分 为 若干 个 相互 重 倒 的 子 阵 列 ， 各 子 阵 列 的 阵列 流 形 相 
同 ， 各 子 阵 列 的 协 方差 矩阵 可 以 相 加 后 进行 平均 运算 ， 取 代 
原来 意义 上 的 协 方差 矩阵 ， 从 而 实现 去 相干 。 
人 


J 


- "| 后 疝 汪 
< ul) [和 
< (0) | 
是 2 m m+1 M-1 M 
#0 
0 
< : "| 
前 向 于 洽 #0 + 


图 2 双向 平滑 算法 子 阵 划分 图 

Fig.2 Bidirectional smoothing algorithm subarray division diagram 
如 图 2 所 示 , 均匀 线 阵 含有 mM 个 阵 元 , 分 成 相互 交错 的 

工 个 子 阵 ， 每 个 子 阵 有 普 个 单元 ， 其 中 四 =M-ZL+1l。 
均匀 线 阵 接收 到 的 信号 矩阵 为 
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下 (四 =[2( 仆 , (1) ,Xm (0 (12) 
根据 图 2 的 分 割 ， 可 得 第 /个 前 向 子 阵 收 到 的 信号 矩阵 
为 
Xi =[% (7), Xn (1), Xi (Dye, Na OF (13) 
Xi =ADTS(1)+N,(t), 1<1<L (14) 
其 中 :DD 是 一 个 对 角 阵 ， D=diag (wv2,…,vp)， Vi=exp(-j@)。 
D=dias(e ,ee,.. ,© | (15) 
则 第 :个 前 向 子 阵 数 据 协 方差 矩阵 为 

RN =E(Xi (1) x7 (CD )=4DDR CD ) 42+o27 (16) 
利用 所 有 的 数据 协 方差 矩阵 求 矩 阵 ， 即 为 前 向 空间 平滑 

协 方差 矩阵 玉 为 
RW -TR (17) 

姑 为 子 阵 阵 元 数目 m>N ， 则 当 L>WN 时 前 向 空间 平滑 数 
据 协 方差 矩阵 R' 是 满 秩 的 N09。 分 析 均 匀 线 阵 的 倒序 阵 ， 同 
理 可 得 后 向 平滑 子 阵 的 接收 数据 为 

XP =[ nm (GEMORE A (18) 
KP =AD"' (DY1S(1)) + NI(t) 1<1<L (19) 
则 每 个 后 向 平滑 子 阵 的 协 方差 矩阵 为 
Rr =E(XP(t) XP (1) )=AD™'R; CD) 42+a27 (20) 
同样 ， 后 向 空间 平滑 的 协 方差 矩阵 为 
及 -Tm CD) 
其 中 : R 是 有 的 共 恩 倒序 阵 。 根 据 共 恩 倒序 不 变性 定义 前 
后 向 平滑 协 方差 矩阵 为 
R=3(R! +R’) (22) 
固 倒 序 不 变性 的 优点 在 于 增加 子 阵 数 目 ， 改 变 了 基于 
天 线 阵 列 协 方差 矩阵 的 特征 分 解 类 DOA 算法 的 局 限 性 。 

下 面 将 前 后 向 数据 协 方差 R 应 用 于 改进 的 IMUSIC 算法 
进行 相干 信 源 DOA 估计 。 改 进 的 IMUSIC 算法 通过 对 接收 
信号 协 方差 矩阵 进行 重 构 ， 使 信号 协 方差 矩阵 的 秩 恢 复 为 
rank(R)=p ， 从 而 能 有 效 地 估计 出 信号 的 DOA。 

令 J 为 MxM 反 单 位 矩阵 ， 即 

0 0 .1 
1 03) 
1 0 0 J 
在 此 令 改 进 后 的 信号 相关 矩阵 为 
R=R+JR'J (24) 
其 中 : R’ 为 及 的 共 恩 ,此 时 玉成 为 Hermite 的 Toeplitz 矩阵 ， 
R 是 RK 的 无 偏 估 计 。 
再 对 R 进行 奇异 值 分 解 得 到 噪声 子 空间 VV: 
[U,S,V]=svd(R) (25) 

取 V=U(:;,p+1:M) 构 造 噪声 子 空间 WV。 用 低 秩 和 矩阵 代 奉 

满 秩 和 矩阵 R ， 令 
S(M-2,M-2)=0,S(M -1,M-1)=0,S(M,M)=0 (26) 


Ss=S,R,=U*Ss*V 


(27) 
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个 秩 比 较 低 的 矩阵 来 近似 
该 问题 可 转换 为 一 个 矩阵 分 


对 于 大 规模 的 矩阵 ， 可 以 用 
己 知 矩阵 ， 即 矩阵 的 低 秩 近似 。 
解 的 问题 ， 即 将 一 个 矩阵 分 解 为 


其 中 :el 为 谱 范 数 或 者 Frobenius 范 数 ; 


(28) 


R=R,+Q 
其 中 : R; 为 低 秩 和 矩阵 ，2 为 稀 疏 和 矩阵 。 算 阵 低 秩 分 解 问题 可 
转换 为 下 面 的 优化 问题 : 
Min am (ER,) = | 及 一 已 外 


(29) 


f 为 矩阵 RE 的 秩 通 


对 新 的 协 方 差 矩阵 RE 进行 奇异 值 分 解 ， 重 复 上 述 步 又 ， 

再 一 次 构造 噪声 子 空间 V。 
[UU,SS,VV]= svd(R,) (30) 
取 V=UU(:,P+1:M) 作为 与 噪声 特征 值 对 应 的 特征 向 量 ， 


即 噪声 子 空间 WV。 


最 后 将 两 次 得 到 噪声 子 空间 矢量 进行 取 平 均值 ， 即 


VU=(V +V,)/2 
前 数 的 空间 谱 ， 即 
Busic (0) = 2 二 2 
[a” (9)VUVU "a(0)] la” (9)vUO 


诗 
奢 
区 | 


Pe (9) 3 


下 : 


模型 由 式 〈1) 所 示 。 


b) 按 照 前 后 向 空间 平滑 技术 对 天 线 阵列 进行 划分 ， 求 的 
前 向 数据 信号 矢量 X/ 和 后 向 数据 信号 矢量 X? 。 根 据 式 (16) 


和 “20) 分 别 求 出 对 应 的 协 方差 矩阵 丸和 展 。 
c) 根 据 步骤 b) 中 得 到 的 数据 矢量 ， 分 别 求 出 
的 数据 协 方差 矩阵 的 平均 值 R* 和 R? 。 


由) 由 Ri 与 R* 的 均值 求 前 后 向 数据 协 方差 矩阵 Rw 。 
6) 令 J 为 MxM 反 单 位 矩阵 ,通过 对 接收 信号 协 方差 矩阵 
使 得 R 成 为 Hermite 的 Toeplitz 


R 进行 重 构 ， 根 据 式 (24) 
矩阵 。 
站 再 对 RR 进行 奇异 值 分 解 得 到 噪声 子 空间 。 
外 然后 用 一 个 低 秩 矩阵 代 莹 满 秩 甜 阵 


FR ， 构造 
协 方差 矩阵 R, 。 重 复 上 述 步骤 ， 并 重 构 噪声 子 空间 WV。 


G32) 


峰值 可 以 通过 改变 9 值 来 进行 入 射 信号 DOA 估计 


根据 以 上 的 讨论 ， 混 合 型 MUSIC 算法 步骤 可 以 总 结 如 


a) M 个 阵 元 组 成 的 天 线 阵列 ， 则 第 1 次 采样 点 数 的 信号 


工 个 子 阵 


一 个 新 的 


通过 计算 噪声 子 空间 的 平均 值 ， 可 得 到 前 后 向 平滑 协 


方差 矩阵 VD 。 


iD 构造 函数 的 空间 谱 Bsc (9) 。 


3 ”仿真 和 分 析 


为 了 更 直观 地 展现 算法 性 能 ， 采 用 MATLAB 对 以 上 算 


ey 
对 


到 的 噪声 为 平稳 零 均 值 高 斯 白 噪 声 , 方差 为 0 且 与 
日 关 。 

实验 1 仿真 模型 为 8 
d=4/2，3 个 信号 源 分 别 以 来 波 方向 -45”、0" 


St 


阵 元 均匀 直线 阵列 ， 阵 元 间 吕 
、30"。 入射 
到 天 线 阵列 ，0"” 和 30” 信号 间 保 持 完 全 相干 的 关系 ， 信 噪 


法 进行 仿真 。 假 设 信 号 满足 远 场 条 件 和 窄带 条 件 ， 阵 元 接收 


目标 源 不 


ial 


比 均 为 10 dB， 采 样 数 为 1 024。 实 验 仿 真 结果 如 图 3 所 示 。 
图 3 明显 地 可 以 得 出 ， 如 果 信 号 是 相干 的 ， 则 经 


MUSIC 算法 的 性 能 严重 下 降 ， 无 法 分 辨 信号 源 的 DOA; 男 


外 三 种 算法 可 以 从 空间 谱 1 


线 上 看 到 三 个 明显 的 峰值 极点 ， 
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并 且 都 对 应 着 精确 的 入 射 角 ， 其 中 混合 型 MUSIC 算法 的 谱 
值 明 显 大 于 其 他 算法 ， 说 明 三 种 算法 均 能 处 理 相 干 信号 但 混 
合 型 MUSIC 算法 性 能 更 优越 。 


”经典 MUSIC 算 法 

160 - -前 后 向 平滑 MUSIC 算 法 | 
9 全 改进 的 IMUSIC 算 法 

性 -混合 型 MUSIC 算 法 | 


-100 -80 -60 -40 -20 20 40 60 80 100 


入 有 角度 
3 _ 信 源 方向 为 [-45%0°30°] 仿 真 结果 
Fig.3 Source direction is [-45° 0° 30°] simulation result 
实验 2 ”仿真 模型 为 8 阵 元 均匀 直线 阵列 ， 阵 元 间距 
4=4/2， 三 个 信号 源 分 别 以 来 波 方向 为 [-2”0”2”] 入 射 到 
天 线 阵列 , 信号 间 保 持 完全 相干 的 关系 , 信 噪 比 均 为 30 dB， 
采样 数 为 1024。 实 验 仿真 结果 如 图 4 所 示 。 


180 1 


吗 


hn 


L | 
*- 经 典 MUSIC 算 法 
-前 后 向 空间 平滑 算法 | 
个- 改进 的 IMUSIC 算 法 || 
信 混 合 型 MUSIC 算 法 


-40 -20 20 40 60 80 100 


图 4 角度 为 [-2”0”2”] 仿 真 结果 
Fig. 4 Angle is [-2° 0° 2°] simulation result 

图 4 结果 可 以 得 出 ， 在 信号 间隔 很 小 且 存 在 相干 关系 
时 ， 经 典 的 MUSIC 算法 和 改进 的 IMUSIC 算法 完全 失效 。 
在 同等 条 件 下 ， 混 合 型 MUSIC 算法 较 之 前 后 向 空间 平滑 算 
法 ,更 能 精确 地 分 辩 相 邻 间隔 的 相干 信 源 。 

实验 3 ”仿真 模型 为 8 阵 元 均匀 直线 阵列 ， 阵 元 间距 
4d=4/2， 三 个 信号 源 分 别 以 来 波 方向 -53” 、0"” 、5” 入 射 到 
天 线 阵 列 ， 信 和 号 间 保 持 完 全 相干 的 关系 ， 同 时 设 定 信 噪 比 
SNR 分 别 在 30dB、20 dB、10 dB、5 dB， 采 样 数 为 1 024 的 
情况 下 对 四 种 算法 进行 仿真 。 实 验 仿真 结果 如 图 5~8 所 示 。 


60r L T ] 
全 30dB 经 典 MUSIC 算 法 
人 -20dB 经 典 MUSIC 算 法 
具 10dB 经 典 MUSIC 算 法 1 
-9-5dB 经 典 MUSIC 算 法 


| 
10 20 30 40 50 


-50 -40 -30 : -10 0 
入 射 角 / 度 


图 5 在 不 同 的 信 噪 比 下 ， 经 典 MUSIC 算法 仿真 结果 
Fig. 3 _ Simulation results of classical MUSIC algorithm at different 


signal-to-noise ratios 
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120 L f | 
\ 全 -30dB 改 进 的 IMUSIC 算 法 
1 时 4 -20dB 改 进 的 IMUSIC 算 法 
100 pd -#-10dB 改 进 的 IMUSIC 算 法 1 
上 | -5dB 改 进 的 MUSIC 算 法 
| 
1 


谱 峰 /dB 


入 射 角 / 度 
图 6 在 不 同 的 信 噪 比 下 ， 改 进 的 IMUSIC 算法 仿真 结果 


Fig. 6. Improved IMUSIC algorithm simulation results at different 


signal-to-noise ratios 


180 | T T T ] 
全 30dB 前 后 向 空间 平滑 算法 
160 上 人 20dB 前 后 向 空间 平滑 算法 1 
宙 一 -10dB 前 后 向 空间 平滑 算法 | 
-5dB 前 后 向 空间 平滑 算法 
120 上 | 
四 1o0 - 
次 
小 80 上 


-50 -40 -30 -20 -0 0 10 20 30 4 50 
入 射 角 / 度 


图 7 在 不 同 的 信 噪 比 下 ， 前 后 向 空间 平滑 算法 仿真 结果 
Fig.7 Simulation results of forward and backward spatial smoothing 


algorithm under different signal-to-noise ratios 


250 T T 1 
全 30dB 混 合 型 MUSIC 算 法 
-20dB 混 合 型 MUSIC 算 法 

| 10dB 没 合 型 MUSIC 算 法 | 


-5dB 混 合 型 MUSIC 算 法 


-50 -40 -30 -20 -10 


0 
入 射 角 / 度 
图 8 在 不 同 的 信 噪 比 下 ， 混 合 型 MUSIC 算法 仿真 结果 
Fig.8 _ Simulation results ofhybrid MUSIC algorithm under different 


signal-to-noise ratios 
图 5~8 结果 可 以 得 出 ， 随 着 信 品 比 的 不 断 递 减 ， 算 法 
的 分 辨 能 力 都 出 现 不 同 程度 的 降低 。 当 SNR=30 时 ， 改 进 的 
IMUSIC 算法 、 前 后 向 空间 平滑 算法 和 混合 型 MUSIC 算法 都 
能 对 信号 实现 DOA 估计 ; 当 SNR=20dB 或 10 dB 时 ， 改 进 
的 IMUSIC 算法 已 经 无 法 分 辨 信号 ， 前 后 向 空间 平滑 算法 出 
现 估计 偏差 , 此 时 的 混合 型 MUSIC 算法 依然 可 以 正确 分 辨 ; 
当 SNR=5dB 时 ， 在 这 种 信号 间隔 很 小 ， 信 噪 比 低 且 信和 号 间 
完全 相干 的 情况 下 ， 混 合 型 MUSIC 算法 还 能 很 好 地 估计 出 


图 9~11 所 示 。 


100 


90 下 


80r 


70r 


60r 


50 


成 功 概率 % 


ChinaXiv 合 作 期 刊 


录用 定稿 张 权 ， 等: 基于 混合 型 MUSIC 算法 对 相干 信 源 DOA 估计 第 37 卷 第 5 期 


句 平滑 算法 而 言 平滑 次 数 为 6)。 信 噪 比 SNR 从 -10 dB ~ 30 
dB 变化 时 ， 得 到 四 种 算法 随 信 噪 比 变 化 的 估计 性 能 曲线 如 


一 米 一 MUSIC 算 法 

一 疡 一 改进 的 IMUSIC 算 法 
一 © 一 前 后 向 空间 平滑 算法 
一 今 一 混合 型 MUSIC 算 法 


对 


9 


Fig. 9 Shows relationship between success probability and 


1 1 1 1 1 
5 10 15 20 25 30 
信 噪 比 〈dB) 


成 功 概 率 与 信 噪 比 的 关系 


signal-to-noise ratios 


一 米 一 MUSIC 算 法 

一 请 一 改进 的 IMUSIC 算 法 
一 合 一 前 后 向 空间 平滑 算 法 
一 售 一 混合 型 MUSIC 算 法 


10 


图 10 


信 噪 比 (dB》 


估计 偏差 与 信 品 比 的 关系 


Fig. 10 Shows relationship between estimated deviation and 


signal-to-noise ratios 


0.5 
一 米 -MUSIC 算法 

0.45 -此 -改进 的 IMUSIC 算 法 
一 @ 一 前 后 向 空间 平滑 算法 

0.4 一 令 一 混合 型 MUSIC 算 法 | ] 


图 11 


15 20 25 30 35 
信 噪 比 (dB) 


估计 方差 与 信 噪 比 的 关系 


Fig. 11 Shows relationship between estimated variance and 


signal-to-noise ratios 


由 仿真 结果 可 以 得 出 ， 图 9 说 明 随 着 信 噪 比 的 增加 ， 四 


估计 概率 最 低 。 图 10 


种 算法 的 估计 概率 都 有 明显 提高 ; 但 总 体 而 言 ， 混 合 型 
MUSIC 算法 在 低 信 噪 比 情况 下 成 功 估计 概率 最 高 , 其 次 是 前 
后 向 空间 平滑 算法 ， 而 经 典 MUSIC 算法 在 低 信 噪 比 情况 下 


中 说 明 在 信 噪 比较 大 时 ,四 种 算法 的 估 


言 号 的 DOA, 实验 证 明了 混合 型 MUSIC 算法 在 面 对 不 同 信 
噪 比 时 具有 较 高 的 性 能 及 稳定 性 。 

实验 4 仿真 模型 为 8 阵 元 均匀 直线 阵列 ， 阵 元 间距 
d=4/2， 两 个 信号 源 分 别 以 来 波 方向 0" 、5” 入 射 到 天 线 阵 
列 ， 信 号 间 保 持 完 全 相干 的 关系 ， 快 拍 数 为 100。 选 取 子 阵 
的 阵 元 数 均 为 6( 即 对 前 向 平滑 算法 而 言 平滑 次 数 为 3, 对 双 


计 偏 差 都 比较 小 ， 即 


较 高 的 估计 精度 ， 而 在 低 信 品 比 情况 


下 ， 混 合 型 MUSIC 算法 依旧 有 较 低 的 估计 偏差 。 图 11 说 明 


混合 型 MUSIC 算法 的 估计 标准 差 最 小 ， 也 就 说 明 混 合 型 
MUSIC 算法 估计 的 稳定 性 较 好 ， 其 他 几 种 方法 相差 不 大 。 


综 上 所 述 ,混合 型 MUSIC 算法 的 估计 信 噪 比 门限 更 低 ， 


估计 偏差 和 估计 标 ; 


差 都 较 低 ， 具 有 较 高 的 估计 精度 和 稳定 


录用 定稿 张 ” 权 ， 


性 ， 特 别 是 当空 间 信 源 方 位 相距 较 近 、 信 噪 比较 低 时 ， 混 合 


型 MUSIC 算法 性 能 的 优越 性 会 更 加 突出 。 


4 ”结束 语 

本 文 提出 了 一 种 用 于 相干 信 源 DOA 估计 的 混合 型 
MUSIC 算法 。 该 算法 实质 上 是 将 前 后 向 空间 平滑 算法 和 改进 
的 IMUSIC 算法 进行 融合 ， 通 过 前 后 向 平滑 技术 对 天 线 阵 列 
进行 预 处理 ， 并 将 得 到 的 新 协 方差 矢量 矩阵 应 用 于 改进 的 
IMUSIC 算法 进行 信号 数据 处 理 分 析 ， 最 后 根据 数值 仿真 分 
析 ， 将 混合 型 MUSIC 算法 与 前 后 向 空间 平滑 算法 和 改进 的 
IMUSIC 算法 进行 对 比 。 仿真 结 果 表 明 , 混合 型 MUSIC 算法 
有 较 高 的 测 向 分 辩 力 ， 能 够 有 效 估计 出 独立 信号 源 和 相干 
言 号 源 的 来 波 方向 ， 并 且 在 信 噪 比较 低 、 信 和 号 间隔 很 小 且 存 
在 强 相 关 时 ， 混 合 型 MUSIC 算法 也 能 很 好 地 估计 出 信号 的 
DOA， 验 证 了 本 文 提 出 的 混合 型 MUSIC 算法 的 高 分 辨 率 和 
高 性 能 。 
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